Communication savings with ghost cell expansion for domain decompositions of finite difference grids
نویسندگان
چکیده
We investigate trade-offs between communications frequency and bandwidth, using the ghost cell expansion method (GCE) with domain decomposition for parallel evaluation of finite difference schemes for PDEs. GCE delays communication of boundary values between local domains for a small number of timesteps. This delay reduces communication pressure and aggregates messages into larger packets, thus optimizing for the high latency, high bandwidth interconnects common in modern parallel computing platforms. However, choosing the optimal ghost cell expansion level requires a careful balance between the cost of local computation per timestep, and network bandwidth (which may limit the effectiveness of increasing message size). Furthermore, for stencils wider than one cell on each side, GCE is an approximation which may impair numerical accuracy. For now, we neglect numerical aspects of GCE. We use a LogP model with bandwidth to show that for the platforms tested, the technique only has potential for reducing the amortized cost of communication in the 1-D and 2-D cases. The model suggests a declining payoff for increasing expansion level in the 2-D case. Experimental results with a test problem verify some of these predictions. 1 Ghost cell expansion Ghost cell expansion, proposed by Ding and He [6], is a form of domain decomposition with overlapping boundaries. Ding and He formulate it for finite difference methods used to solve PDEs with a timestep iteration. It accomodates either time-dependent PDEs (parabolic or hyperbolic) or elliptic problems E-Mail: [email protected] E-Mail: [email protected]
منابع مشابه
Water hammer simulation by explicit central finite difference methods in staggered grids
Four explicit finite difference schemes, including Lax-Friedrichs, Nessyahu-Tadmor, Lax-Wendroff and Lax-Wendroff with a nonlinear filter are applied to solve water hammer equations. The schemes solve the equations in a reservoir-pipe-valve with an instantaneous and gradual closure of the valve boundary. The computational results are compared with those of the method of characteristics (MOC), a...
متن کاملA Supra-Convergent Finite Difference Scheme for the Poisson and Heat Equations on Irregular Domains and Non-Graded Adaptive Cartesian Grids
We present finite difference schemes for solving the variable coefficient Poisson and heat equations on irregular domains with Dirichlet boundary conditions. The computational domain is discretized with non-graded Cartesian grids, i.e., grids for which the difference in size between two adjacent cells is not constrained. Refinement criteria is based on proximity to the irregular interface such ...
متن کاملA Novel Structure for Optical Channel Drop Filter using Two-Dimensional Photonic Crystals with Square Lattice
In the present paper a novel structure for optical channel drop filter (CDF) based on photonic crystal ring resonator with circular core has been proposed. In order to design the proposed CDF, the plan wave expansion (PWE) method is applied for calculation of band structure and photonic band gap while the transmission characteristics of proposed CDF have been calculated using the finite differe...
متن کاملFinite-Difference Time-Domain Simulation of Light Propagation in 2D Periodic and Quasi-Periodic Photonic Structures
Ultra-short pulse is a promising technology for achieving ultra-high data rate transmission which is required to follow the increased demand of data transport over an optical communication system. Therefore, the propagation of such type of pulses and the effects that it may suffer during its transmission through an optical waveguide has received a great deal of attention in the recent years. We...
متن کاملMimetic Finite Difference Modeling Of 2D Elastic P-SV Wave Propagation
Recently, efforts have been made by Kristek et al.(2002) toward the implementation of a zero-traction boundary condition for an elastic medium by using high-order staggered-grid finite-difference modeling and avoiding symmetry conditions, vacuum formulations, or, other approaches that require grid points located beyond the physical boundary (ghost points). In this work, a new set of numerical d...
متن کامل